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ABSTRACT 


Late type giants (RGB/AGB stars) will be important 
tracers of the Galactic morphology and evolution in the 
framework of Gaia, as they are intrinsically bright and 
thus can probe distant stellar populations or those ob¬ 
scured by interstellar extinction. A realistic represen¬ 
tation of their atmospheres and spectra with stellar at¬ 
mosphere models is thus of crucial importance, both 
for the design and optimization of Gaia instruments, 
as well as the interpretation of provided astrophysical 
data. Our analysis of synthetic photometric colors of 
late-type giants based on PHOENIX, MARCS and ATLAS 
model atmospheres indicates a general agreement be¬ 
tween the current theoretical predictions and observations 
in the framework of stationary 1-D model atmospheres. 
Presently available models allow temperature determina¬ 
tions of RGB/AGB stars to an accuracy of ~ ±100 A. In 
an exploratory study we try to quantify possible resid¬ 
ual systematic effects due to the approximations made 
in 1-D models using full 3-D hydrodynamical models. 
We find that differences in broad-band photometric col¬ 
ors calculated with 1-D and 3-D models are significant, 
translating to the offsets in effective temperature of up to 
A T e ff ~ 70 K. Clearly, full 3-D hydrodynamical models 
will help to alleviate such ambiguities in current theoret¬ 
ical modeling. Additionally, they will allow to study new 
phenomena, to open qualitatively new windows for stellar 
astrophysics in the Gaia-era. 
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1. IMPORTANCE OF LATE-TYPE GIANTS IN 
THE CONTEXT OF GAIA 


Gaia will be an excellent instrument for studying Galac¬ 
tic stellar populations, especially in terms of their detailed 
formation histories, as it will provide parallaxes, proper 
motions, radial velocities, and astrophysical parameters 


(effective temperatures, metalicities, gravities and red¬ 
denings) for stars down to Gaia magnitude G = 20, 
which translates to Johnson V ~ 21 for V — I = 2.0 
(or T eS ~ 3700 K) and V ~ 20.5 for V - I = 1.0 
(T e ff ~ 4800 K). This unique combination of kinemat- 
ical and astrophysical information which is not accessi¬ 
ble otherwise will provide an opportunity to investigate 
formation and evolution of the Galaxy in unprecedented 
detail. 

Currently stellar populations are predominantly dated us¬ 
ing the main sequence turn-off point (MSTO) stars which 
is the most reliable method among those available to¬ 
day. This requires precise photometry of typically faint 
main sequence stars, i.e. those located at least one mag¬ 
nitude below the MSTO (at 15 Gyr and solar metalicity 
the MSTO is at My ~ ±4.6). Moreover, photometric 
data should be obtained for a statistically sufficient num¬ 
ber of stars in a given target population. Surprisingly, 
these two requirements appear to be rather demanding in 
the context of Gaia, since already a mild extinction in 
the galactic disk of e.g. Ay ~ 0.7mag/kpc will make 
MSTO stars too faint to be detected with Gaia at rather 
short distances. The limiting distance to which stellar 
population can be dated with Gaia using MSTO stars will 
vary depending on the number of MSTO stars, interstel¬ 
lar extinction, and other factors, but it will be generally 
confined to 2.0 — 2.5 kpc ( :Bertellil2002h . Clearly, unless 
other suitable and reliable tracers can be identified, the 
capability of Gaia for dating stellar populations will be 
restricted to a few kiloparsec within the solar neighbor¬ 
hood. Late-type giants on the RGB and early-AGB can 
be very useful in overcoming this limitation. Being sig¬ 
nificantly brighter than MSTO stars they can be easily ob¬ 
served to large distances (at solar metalicity the absolute 
magnitude of a typical RGB star located ~ 1.5 mag be¬ 
low the RGB tip is My ~ ±0.5). Hence, late-type giants 
will be extensively employed in the Gaia data analysis as 
they can provide information about distant stellar popula¬ 
tions not accessible by using fainter main sequence stars. 

Until now, however, late-type giants have been scarcely 
used for dating stellar populations. Firstly, the obser- 
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vational scatter in their photometric colors is typically 
large compared to the spread of isochrones correspond¬ 
ing to different ages on the RGB/AGB, thus uncertainties 
in derived ages are large. Secondly, the evolution on the 
RGB/AGB is still rather uncertain, as theoretical models 
are dependent on a number of poorly constrained model¬ 
ing parameters such as the mixing-length parameter and 
mass loss rate. 

However, for Gaia there is evidence that effective temper¬ 
ature determinations of late-type giants may be feasible at 
an accuracy level of 1 — 2%, mostly from medium-band 
photometry (see Jordi & H0g, this volume). The quoted 
error level corresponds to the formal error resulting from 
the fitting procedure with theoretical model atmospheres. 
At least for some stars the medium-band photometry will 
be supplemented by ground-based data obtained by mas¬ 
sive future spectroscopic (e.g. RAVE) and photometric 
(e.g. VST, VISTA, UKIDSS) surveys. Near-infrared pho¬ 
tometry will be especially important in this context, as 
it will allow to put even tight er constrains on the effec¬ 
tive temperatures (see iKucinskas et al.ll2003l K03) and 
interstellar extinction. This will allow to obtain suffi¬ 
ciently precise astrophysical parameters of RGB/AGB 
stars (cr(Teff) ~ 50 — 100 K - see Sect. 2), which, in 
turn, may provide age estimates at an (formal) accuracy 
level of ~ 20% or better. 

The most essential ingredient missing here today is the 
availability of reliable evolutionary models, which would 
provide a solid basis for precise and unambiguous dat¬ 
ing of stellar populations with RGB/AGB stars. At least 
partly this is due to the inadequate treatment of multi¬ 
dimensional phenomena (especially convection) within 
the classical 1-D models. It is conceivable that at least 
some of these problems will be solved with full 3-D hy- 
drodynamical models, which may be available for routine 
use within ~10 years from now. Being constructed to ac¬ 
count for time-dependent and non-spherical phenomena, 
3-D hydrodynamical models may help to alleviate the 
majority of ambiguities in theoretical modeling of stel¬ 
lar evolution on the RGB/AGB (especially those related 
with convection, mass-loss, pulsations, etc.). 


2. CURRENT STATUS IN THEORETICAL 
MODELING OF LATE-TYPE GIANT AT¬ 
MOSPHERES: SYNTHETIC BROAD-BAND 
PHOTOMETRIC COLORS 


Clearly, a number of issues related to stellar evolution on 
the RGB/AGB have to be clarified before late-type giants 
can be widely used for reliable dating of stellar popula¬ 
tions. Can today’s theoretical model atmospheres provide 
a sufficient backup to obtain precise astrophysical param¬ 
eters of individual giants from observations, in particular 
broad-band photometry? 

To address this question we have performed an extensive 
comparison of synthetic broad-band photometric colors 
of late-type giants with observations, covering a wide 
range of effective temperatures, gravities and metalici- 
ties (T eff = 3000 - 4500 K, log g = 0.5 - 2.5 and 
[M/H] = 0,—1.0,—2.0). Synthetic spectra used for 


this purpose were produced with the most recent ver¬ 
sions of PHOENIX, MARCS and ATLAS model atmo¬ 
spheres (for details see: PHOENIX: Brott & Hauschildt, 
this volume; MARCS: Plez 2003; ATLAS: Castelli & Ku- 
rucz, 2003). Synthetic colors were calculated in the 
Johnson-Cousins-Glass system, using filter definitions 
from lBes selll ill 99 (1) for th e John son-Cousins BVI bands 
and from IBesseil & Breti d 1988b for Johnson-Glass JK 
bands. A carefully selected sample of local late-type 
giants with T e s < 5000 K was used for this compari¬ 
son, with precise determinations of effective temperature 
available from interferometry (errors in T e g typically be¬ 
low 5%). We also made sure that stars in this sample are 
non-peculiar (chemically or otherwise) and non-variable. 
Comparisons were made in T e ff — color and color—color 
planes, with synthetic photome tric colors plotted accord - 
ing to a T e ff — log g relation of iHoudashelt et aft (200(1). 
We included several commonly used T e ff — color rela¬ 
tions in this analysis too (Fig.l; A99: Alonso et al. 1999; 
BaSel2.2: Lejeune et al. 1998; SF00: Sekiguchi & 
Fukugita 2000; VC03: Vandenberg & Clem 2003; H00: 
Houdashelt et al. 2000). Here we present some results 
for solar metalicity; a full discussion covering also sub¬ 
solar metalicities will be given elsewhere (Kucinskas et 
al. 2005, A&A, in preparation). 

Figure 1 shows that there is a general agreement between 
the synthetic broad-band colors predicted by today’s sta¬ 
tionary 1-D model atmospheres and the observations, es¬ 
pecially in the T e ff — (V — K) and T e ff — (J — K) planes. 
Moreover, there is a remarkable agreement among the 
synthetic colors calculated with different stellar atmo¬ 
sphere models, translating typically into an uncertainty 
in T e ff of ±50 K over a large range of effective tem¬ 
peratures, gravities and metalicities. This indicates that 
presently available models may allow to obtain effec¬ 
tive temperatures of late-type giants with a precision of 
~ ±100 K, using a single broad-band photometric in¬ 
dex like V — /, or V — K. The accuracy can still be 
improved if optical and near-infrared broad-band fluxes 
(such as Johnson’s BVRIJHK) are fitted simultaneously 
with synthetic fluxes to obtain an estimate of T e g (K03). 
Typically, observed fluxes at 5-6 wavelength points are 
sufficient to achieve a formal fitting accuracy (standard 
deviation) of cr(T e ff) ~ 50 K. This considerably reduces 
the observational scatter on the RGB/AGB, allowing to 
obtain ages with a precision of at least ±20% (see K03 
for an example with LMC/SMS clusters). 


3. TOWARDS HIGHER ACCURACY: FULL 
3-D HYDRODYNAMICAL MODEL ATMO¬ 
SPHERES 


We conducted an exploratory study (Ludwig et al. 2005, 
in preparation) to estimate systematic effects on broad¬ 
band colors related to changes of the mean temperature 
profile and the presence of temperature inhomogeneities 
in a fully time-dependent 3-D hydrodynamical model at¬ 
mosphere relative to a 1-D standard one. We consid¬ 
ered a prototypical late-type giant with T e ff = 3700 K, 
log g = 1.0, [M/H\ = 0. It is an interesting example 
since standard models predict that convection is restricted 
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Figure 1 . Empirical and theoretical T e g — color relations for late-type giants in different T e ff — color planes (a-d, top 
panels). Filled circles are non-variable late-type giants, several semiregular variables are highlighted as open circles. 
Thick solid line is a best-fit to the obsen’ed colors in a given T e g — color plane. Several existing T e g — color relations 
are shown as well, together with T e g — color scales constructed using synthetic colors o/PHOENIX, MARCS and ATLAS 
(see text for details). Bottom panels in each figure show the difference between various T e g — color relations and the best 
fit to the obsen’ed data in a given T e g — color plane ( A T e g = T°^ her — T e l g stfit J. 
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Figure 2. Snapshot of the emergent white light inten¬ 
sity during the temporal evolution of convective flows on 
the surface of a red giant (T e g = 3700&, logg = 1.0, 
[M/H] = 0). A typical granular pattern is discernible. 


to the optically thick layers, hence has no influence on the 
atmospheric structure and consequently on the spectro¬ 
scopic or photometric stellar properties. In contrast (see 
Fig. our hydrodynamical model predicts that convec¬ 
tion affects the atmospheric layers, and a typical granu¬ 
lation pattern is present. As it turns out it is the result of 
intense convective overshooting. 

1-D model atmosphere used in this comparison employed 
the same physical input data (opacities, equation of state, 
description of radiative transfer) as the hydrodynamical 
model. The convective energy transport was treated with 
mixing-length theory, turbulent pressure was neglected. 
Both aspects are of minor importance here since neither 
process influences the atmospheric structure of the 1-D 
model. For the 3-D model, direct calculation of broad¬ 
band colors is computationally too demanding, because 
of the complex geometry and time-dependence of the 
convective flow, and a large number of wavelength points 
necessary for a realistic description of the stellar spec¬ 
trum. We instead approximated the hydrodynamical flow 
structure by classifying horizontal locations according to 
their white light emergent intensity, by sorting darker and 
increasingly brighter areas into nine intensity groups. For 
each group we calculated the average (in space and time) 
vertical thermal structure. The resulting nine thermal 
structures were treated as standard plane-parallel model 
atmospheres (the so-called 1.5-D approximation) in the 
subsequent spectral synthesis calculations. The resulting 
radiation fields were added according to the surface area 
occupied by their intensity group. Figure 0 shows the 
magnitude differences between broad-band fluxes of the 
3-D and 1-D model. We find color biases corresponding 
to temperature offsets of up to AT e g ~ 70 K. This is a 
preliminary result obtained for a particular case but it il¬ 
lustrates limitations in the accuracy which we can in prin¬ 
ciple expect when interpreting observations in the frame¬ 
work of hydrostatic 1-D models. By the time Gaia is fly¬ 
ing we expect that hydrodynamical model atmospheres 
are generally available for overcoming such limitations. 

3-D hydrodynamical model atmospheres will also allow 
to account for a number of qualitatively new effects. We 
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Figure 3. Influence of surface granulation on the broad¬ 
band photometric colors of a red-giant, as reflected by 
differences between the predictions of 3-D full hydrody¬ 
namical and classical 1-D model atmospheres. Note that 
differences in color indices - in this representation di¬ 
rectly given by the magnitude differences of the two se¬ 
lected bands - are significant, e.g. A(V — K) ~ 0.2 
would correspond to A T e g ~ 70 K. 


find, for instance, that surface granulation due to convec¬ 
tive overshoot in the upper photospheric layers will in¬ 
duce a luminosity variability on a peak-to-peak level of 
~ 0.01 mag, on a time scales of one week. Addition¬ 
ally, intensity fluctuations on the stellar surface will also 
cause a stochastic movement of the stellar photocenter 
on the same time scale (iSvensson & Ludwial2005t) . Al¬ 
most certainly, such phenomena will be routinely studied 
in 10-15 years from now, thus providing new windows 
for looking inside stars. 
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